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Abstract 

We investigate the ordering of voter model on fractal lattices: Sierpinski Carpets 
and Sierpinski Gasket. We obtain a power law ordering, similar to the behavior of 
one-dimensional system, regardless of fractal ramification. 



1 Introduction 



The Ising model is a well known dynamical model that was investigated in 
complex networks and fractal structures [1,2,3,4,5]. However, aside from that 
model, there are many other possible dynamics, sharing little in common with 
behavior of the Ising model. The voter model is an example of such a model, 
that exhibits different qualities at a very basic level. Unlike the Ising model, 
the voter model has no surface tension and defines a broad universality class 
[6]. While the Ising model dynamics has been studied on fractal lattices [3,4,5] 
little is known about the behavior of voter model in such geometries. 
We have investigated the behavior of the voter model on Sierpinski carpets 
and on Sierpinski gasket. It is known [7] that the evolution of the voter model 
depends on the dimensionality of the lattice. For a large time t the ordering 
process obeys the following equations 



p(t) ~ < 



t~ a , D < 2 

(\nt)-\D = 2 (1) 
1, D > 2 
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Fig. 1. The construction of Sierpinski Carpet (SC) fractal network. Two different 
basic patterns are on the left, and corresponding fractal networks of level 3 and 2 
are on the right. The hatched area of patterns are full positions. 

where p is a fraction of links that form interfaces, i.e. they connect opposite 
spins, D is dimensionality and t is time. The predicted exponent a — 1 — D/2 
was calculated for a lattice with integer dimensionality (-0 = 1). 
We will focus on the question, whether the dynamics on the fractal lattices is 
the same, or are there other rules governing them. 



2 Models 



The voter model is a very simple model of opinion formation. Nodes in the 
network are agents, each one having an opinion. There are only two possible 
opinions, and typically they are considered as +1 and -1, just as Ising spins. 
The dynamic rule is simple — the node opinion changes to an opinion of one 
randomly chosen neighbor. 

The implementation is following: we choose one node at random, and then 
one of its neighbors randomly. The first node assumes the state of the second. 
One time step of the dynamics corresponds to the number of individual node 
updates equal to the number of nodes in the network, so on average each node 
is updated once every time step. 

We investigate the voter model behavior on two fractal networks: Sierpinski 
Carpet (SC) and Sierpinski Gasket (SG). The SC is constructed according to 
a chosen basic pattern. The pattern is a square, divided into nxn squares that 
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Fig. 2. The construction of Sierpinski Gasket (SC) fractal network of level 5. 



could be full or empty (Fig.l). First, single nodes are taken, and arranged into 
the pattern, putting nodes into full positions and skipping empty positions. In 
the next step, the resulting structures are arranged into the same pattern. All 
neighboring nodes in the resulting pattern are connected creating the fractal 
network. The fractal dimension of SC depends on the basic pattern. Classical 
SC has 3x3 pattern with all the squares full except the central one. Such SC 
has a fractal dimension d = |^| « 1.893. We have investigated SC of dimen- 
sions ranging from m 1.6025 to Since it is impossible to numerically 
investigate true, infinite fractals, we will call the number of steps in what the 
network was made a level of fractal. We have investigated SC with different 
patterns, that all have full squares along the pattern edges and empty interior, 
but differ in size. All SC fractals have infinite ramification. 
The ramification is the minimal number of links that one has to remove to sep- 
arate a part of any chosen size from a network. The finite ramification means 
that the structure has some "weak points" where only a finite number of links 
connect together two parts of an infinite network. The infinite ramification 
means that infinite parts of infinite network are connected by infinite num- 
ber of links. For example, a regular square lattice has an infinite ramification, 
while a tree has a finite ramification order. 

The SG network is created in the following way (Fig. 2). Three nodes are taken 
and connected into a triangle. In the middle of each edge a node is created and 
the three new nodes are connected between themselves. This way the whole 
triangle is divided into four smaller ones. In the next step all three non-central 
triangles are treated in the same way, adding nodes in middle of the triangle 
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Interface density evolution for 3 different simulations on SG network 
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Fig. 3. The ordering process in the SG models for three different simulations. The 
data are for fractals of level 8. 

edges and linking them toghether. SG has the fractal dimension |^| 1.5850 
and it possesses a finite ramification. 

While in the case of SC, it was easy to create a general class of SC fractals 
with different fractal dimensions, we are not aware of any generalization of SC 
model that allows easy tuning of fractal dimensions. 



3 Results 

We have investigated ordering of the voter model in SC and SG fractals. 
To measure the disorder, we have used the fraction of interfaces p = I/E, 
where / is a number of interfaces - links connecting nodes with different spins, 
E = N (k) /2 is the total number of links in the network. 

The system orders (Fig. 3) with the interface fraction p decreasing as a power 
of time t. However, due to the finite system size, there are fluctuations around 
the power-law. Since the power-law decay becomes slower with time, the fluc- 
tuations become more significant, and they push the system into a completely 
ordered absorbing state after some time. 

To extract the power-law trend, we have averaged the results of many simu- 
lations, but to avoid the exponential decay due to complete ordering of the 
individual simulations, in a given time step we have averaged only over the 
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Averaged interface density evolution for different networks 
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Fig. 4. The ordering process after averaging over simulations that did not order 
completely. The x marks are for SG, while other symbols are for SC. All the data 
are averaged over 100 network simulations and log-binned. The exponents a exp are 
obtained from the slopes. The SG simulation data has been taken with 10 time steps 
intervals, thus the graph begins later than SC ones. We cannot explain the plateau 
for SG data, but in another simulations with higher level, but lower statistic, the 
plateau was absent. 

simulations that were not completely ordered at that time. This way we have 
circumvented the fluctuations ordering the system and have obtained an ap- 
proximation of an infinite network (Fig. 4). 

We observe the evolution of the interface fraction p in time for networks with 
various fractal dimensions between 1 and 2. We have measured the exponent 
a of that power-law, and compare to the theoretical value [7] (Eq.l). 



4 Conclusions 

The results (Fig. 5) we have obtained show that the theoretical predictions [7] 
(Eq.l) correctly describe the power-law ordering for a voter model on fractals, 
but the exponents obtained through simulations differ significantly from those 
obtained from the analytic formula. 

Moreover, the fact that the voter model in SC and SG behaves in the same 
fashion suggests that the ramification of the fractal does not influence the 
voter dynamics, unlike the Ising model case [3]. 

This work was partially supported by a EU Grant Measuring and Modelling 
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Table 1 

Theoretical and experimental exponents a (see Eq.l) for ordering processes in var- 
ious networks. The results are averaged over 100 individual simulations. The levels 
of the fractals were maximized while keeping a number of nodes that allowed the ac- 
tual simulations to be completed in reasonable amount of time. The regular network 
was a simple 1-dimensional chain with (k) = 4 and periodic boundary conditions. 



Ordering exponent as function of network dimension 




dimension 

Fig. 5. Theoretical and experimental exponents for the ordering processes in various 
networks. The line is the analytic formula (Eq.l), the squares are exponents for 
SC networks, the triangle is exponent for SG network, the circle in the left upper 
corner is the exponent for the one-dimensional network. Error bars are no larger 
than symbol sizes. 
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